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Abstract. We report R C I C light curves of 2 novae in the M31 galaxy which were detected in the four year Nainital 
Microlensing Survey. One of these novae has been tracked from the initial increase in flux while other has been 
observed during its descending phase of brightness. The photometry of the first nova during the outburst phase 
suggests its peak R c magnitude to be about 17.2 mag with a flux decline rate of 0.11 mag day -1 which indicates 
that it was a fast nova. A month after its outburst, it shows reddening followed by a plateau in I c flux. The second 
nova exhibits a bump in R c and J c , possibly about three weeks after the outburst. 
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1. Introduction 

Cataclysmic variables (CVs) are close binary systems con- 
sisting of a white dwarf primary and a late-type main 
sequence secondary star. Novae are a sub-class of cata- 
clysmic variables characterized by the presence of a sud- 
den increase of brightness, called outbursts, due to ther- 
monuclear runway in the envelope of the primary, causing 
the system brightness to increase typically by 10-20 mag. 
These are bright objects which reach up to My ~ -9.0 mag 
at maximum and their rate of decline is tightly correlated 
with their absolute magnitude at maximum (McLaughlin 
1945). The study of novae in external galaxies is impor- 
tant to infer their distances as these objects are one of 
the brightest standard candles up to the Virgo cluster (cf. 
Jacoby et al. 1992 for a review) and tracers of differences 
in the stellar content among galaxies (cf. Van den Bergh 
1988 for a review). 

M31, our nearest large galaxy, has been a target of 
searches for novae since the pioneering work of Hubble 
(1929). Later Arp (1956), Rosino (1964, 1973), Rosino 
et al. (1989), Ciardullo et al. (1987), Sharov & Alksnis 
(1991), Tomaney & Shatter (1992), Rector et al. (1999) 
and Shatter & Irby (2001) have extended the system- 
atic search for novae in M31. In collaboration with the 
AGAPE (Andromeda Gravitational Amplification Pixel 
Experiment) group, we started Cousins R and / pho- 
tometric observations of M31 in 1998 to search for mi- 
crolensing events. Based on the 4 year observations, we 



Send offprint requests to: Y. C. Joshi, 
e-mail: yogesh@upso . ernet . in 



have already reported the discovery of new Cepheids and 
other variable stars (Joshi et al. 2003a). As a microlensing 
survey program is ideally suited to monitor the flux and 
temperature variation during transient events, we have ex- 
tended our search to detect nova outbursts. Here we report 
photometric light curves of two novae detected in the tar- 
get field, one each in 2000 and 2001 observing seasons. 
We show that an increase in flux at longer wavelengths 
a few weeks after the initial rapid decline is a common 
phenomenon in the novae light curves. 

2. Observations 

We have undertaken a program called the Nainital 
Microlensing Survey to detect microlensing events in 
the direction of M31, at the State Observatory, Nainital, 
India, since 1998. Cousins R and I broad band CCD ob- 
servations of M31 were carried out for an ~ 6' x 6' field in 
1998 and an ~ 13' x 13' field during 1999 to 2001 observing 
seasons using the 104-cm Sampurnanand Telescope at the 
f/13 Cassegrain focus. The total integrated observing time 
devoted to the survey ranges from ~ 30 minutes to 2 hours 
each night. The 13' x 13' target field (a 2 ooo = 0' l 43 m 38 s 
and £2000 = +41°09'.l) is centered at a distance of ~ 15 
arcmin away from the center of M31. The average seeing 
during the 141 observed nights spanning 4 years was ~ 2 
arcsec. An overview of the observational detail has been 
given by Joshi et al. (2003a). 

A large database collected during the observing period 
was planned as: 

(a) To search for microlensing events. 
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Fig. 1. A 2 arcmin wide subset of two different R band images taken on January f f , 2000 (left panel) and October 
20, 2000 (right panel). North is at the top and East is to the left. The left frame shows no star at the position marked 
by a circle while right frame shows a star (nova NMS-1) of brightness R ~ 17.2 mag at that position. 



(b) To search for variable stars, particularly Cepheids. 

(c) To search for other transient events e.g. novae, 
where a decrease in brightness of 2 to 5 mag within the 
first two months of the nova evolution can easily be de- 
tected in our observations. 

We have already published a catalog of variable stars 
which contains data for 26 Cepheids and 333 red vari- 
able stars (Joshi et al. 2003a). We are in the process of 
analysing some of the possible microlcnsing candidates 
found in our data set (Joshi et al. 2003b, 2003c). Here 
we report the photometry of two novae detected in the 
target field. 



4. Identification of the novae 

We identified two novae in the target field while searching 
for microlcnsing events using the pixel technique which 
was initially proposed to detect microlcnsing events by 
Baillon et al. (1993). The implementation of the pixel 
technique in our data is described in detail by Joshi et 
al. (2001, 2003c). The two novae detected in our survey, 
one in 2000 and the other in 2001 observing seasons, are 
named nova NMS-1 and nova NMS-2 where NMS is an 
acronym for our project 'Nainital Microlensing Survey'. 
The two novae are individually discussed in the following 
subsections. 



3. Data reduction and photometry 

Standard techniques were used for data reduction using 
MIDAS and IRAF software. The dark current correction 
was not applied due to its negligible contribution dur- 
ing the maximum exposure time of a frame. Cosmic rays 
(CRs) were removed from each frame independently. We 
added all the frames of a particular filter taken on a single 
night and made one frame per filter per night to increase 
the signal to noise ratio. The CCD frames were processed 
using the DAOPHOT photometry routine (Stetson 1987). 
To obtain R and / standard magnitudes, the photomet- 
ric calibration was done using Landolt's (1992) standard 
field SA98, on a good photometric night of 25/26 October, 
2000. A total of 13 secondary stars having 0.09 < (R-I) < 
1.0 were observed over a wide range of airmasses. The typ- 
ical error in magnitudes at ~20.0 mag level is ~ 0.10 and 
0.15 mag in R and / bands respectively. In order to use ob- 
servations in non-photometric conditions, differential pho- 
tometry was done assuming that the errors introduced due 
to colour difference between nova and comparison stars 
were much smaller than the zero point errors. Further de- 
tails of the photometric calibration has been given in Joshi 
et al. (2003a). 



4.1. Nova NMS-1 

The nova NMS-1 having celestial coordinates a20oo = 
h 42 m 57M and <5 2 ooo = +41°07'15".7 was reported in 
the IAU circular by Donato et al. (2001). When we started 
observations at 18:54 UT on October 18, 2000, the nova 
was still brightening. We followed it through the 2000- 
2001 observing season. Fig. 1 shows two images where the 
nova is unresolved and at maximum brightness. 

The photometry of nova NMS-1 has been carried out 
using the PSF profile fitting technique. Since it was suf- 
ficiently bright during the early phase of its eruption and 
gradually increased its brightness, we carried out photom- 
etry on each individual frame of the observations of 18 and 
20 October, 2000 (see Table 1). The light curve of nova 
NMS-1 is shown in Fig. 2. The brightness of the nova 
increased by ~ 0.10 mag in the R band from 18:54 UT 
to 19:17 UT and by ~ 0.23 mag in / band from 19:35 
UT to 20:07 UT. From a comparison of individual frames 
taken on that day, it appears that we have captured nova 
NMS-1 during its peak brightness. The total brightness of 
the nova increased by ~ 1.5 mag during October 18-20, 
2000 followed by an almost exponential decay. Its rate of 
decline, Vd, defined as the rate in mag day -1 at which a 



Y. C. Joshi et al.: Photometric study of two novae in M31 



3 



Table 1. The R, I magnitudes of the nova NMS-1 during outburst. The fractional value for the Julian date is given for those 
points where we calibrated the magnitude for each individual frame of that night. 



J.D. 


R 


J.D. 


/ 


J.D. 


R 




(+2450000) 


(mag) 


(+2450000) 


(mag) 


(+2450000) 


(mag) 


(mag) 


1836.287 


18.69±0.02 


1836.316 


18.38+0.03 


1869 


19.64+0.03 


19.32+0.06 


1836.292 


18.64±0.02 


1836.321 


18.38+0.02 


1870 


19.67+0.02 


19.23+0.05 


1836.296 


18.60±0.02 


1836.327 


18.36+0.02 


1872 


19.69+0.03 


19.39+0.05 


1836.299 


18.57±0.02 


1836.333 


18.36+0.02 


1874 


19.80+0.02 


19.42+0.05 


1836.303 


18.59±0.02 


1836.338 


18.15+0.05 


1877 


19.84+0.03 


19.34+0.07 


1838.213 


17.19±0.01 


1838.236 


16.97+0.01 


1878 


19.93+0.02 


19.24+0.06 


1838.218 


17.16±0.01 


1838.245 


16.95+0.01 


1879 


20.00+0.03 


19.44+0.08 


1838.224 


17.19±0.01 


1838.250 


16.94+0.01 


1881 


19.97+0.02 


19.70+0.06 


- 


- 


1838.255 


16.96+0.01 


1884 


20.13+0.02 


19.55+0.05 


- 


- 


1838.259 


16.94+0.01 


1888 


19.97+0.07 


19.35+0.08 


1843 


17.74±0.01 




17.27+0.01 


1890 


20.18+0.06 


19.37+0.05 


1845 


18.09±0.02 




17.67+0.02 


1905 


20.88+0.06 


19.72+0.05 


1851 


18.69±0.02 




18.44+0.03 


1915 


20.60+0.12 


19.40+0.12 


1853 


18.84±0.02 




18.54+0.05 


1918 




19.12+0.08 


1855 


19.04±0.02 




18.71+0.04 


1921 


20.82+0.08 


19.40+0.06 


1862 


19.20±0.03 




18.89+0.04 


1924 


20.96+0.07 


19.56+0.04 


1863 


19.35±0.03 




19.03+0.04 


1925 


20.81+0.06 


19.52+0.04 


1866 


19.57±0.02 




19.29+0.07 


1926 


20.85+0.08 




1867 


19.61±0.03 




19.09+0.04 


1927 


21.28+0.09 


19.50+0.07 


1868 


19.68±0.02 




19.19 + 0.06 


1929 


21.63+0.11 
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24.49 mag for M31 and a total extinction of 0.63 mag and 
0.47 mag towards our observed direction in the R and / 
bands respectively (cf. Joshi et al. 2003a) . The correlation 
between observed R ma x and rate of decline is consistent 
with that illustrated by Capaccioli et al. (1989) for novae 
in M31. 
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Fig. 2. R, I light curves of the nova NMS-1 during out- 
burst where filled and open circles represent R and I mag- 
nitudes respectively. The arrow indicates the time where 
the nova falls 2 mag below its maximum brightness. 

nova drops to two magnitudes below maximum brightness, 
is estimated to be: 



Vd,R ~ 0.11 mag day 1 
Vd,i ~ 0.11 mag day -1 



(1) 
(2) 

for R and I bands respectively. The ob- 
served values of Vd suggest that the nova 
NMS-1 was a fast nova. The maximum magnitude 
of the nova NMS-1 at peak brightness is Mr (max) 
~ —7.96 and Mj (max) ~ —8.02 mag in R and / bands 
respectively. Here we consider a true distance modulus of 



4.2. Nova NMS-2 

The nova NMS-2 was already reported in the IAU circular 
by Li (2001). Fig. 3 shows two images where the nova is at 
unresolved and maximum brightness phase. The celestial 
coordinates of the nova NMS-2 are a 2 ooo = 00' l 43 m 03 s .3 
and #2000 = +41°12'10 .8. Since this nova was situated 
at the edge of our target field, it could not be observed 
in all the images. When we started our observation on 
October 12, 2001, it was already in a descending phase of 
brightness. 

Photometry of the nova was carried out during the 
brighter phase of outburst. We give its R and / magni- 
tudes in Table 2 and the light curve is shown in Fig. 4. A 
variation of more than 1.5 mag in R band is seen during 
the first 35 days of its observations. For the nova NMS- 
2, our photometry is insufficient to characterize its speed 
class since we do not have any information of its maximum 
brightness. The brightest magnitude of the nova NMS-2 
in our observations is estimated to be ~ 17.7 mag in both 
R and / bands indicating that it must be a star brighter 
than Mr, Mr ~ —7.4 mag during its peak brightness. 

The decay rate for the nova NMS-2 during our obser- 
vation period is estimated to be: 

dM(R) 



0.03 ± 0.003 mag day" 



(3) 



dt 
dM(I) 
df~ 

for R and / bands respectively. The small value of the 
decay rate for the nova NMS-2 suggests that either it was 



0.05 ± 0.004 mag day" 



(4) 
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Fig. 3. A 2 arcmin wide subset of two different R band images taken on January 17, 2001 (left panel) and October 
19, 2001 (right panel). North is at the top and East is to the left. The left frame shows no star at the position marked 
by a circle while right frame shows a star (nova NMS-2) at that position of brightness R ~ 17.7 mag. 



a slow nova or, most likely, we observed it very late after 
it reached its peak brightness. 

A re-brightening in the light curve of the nova NMS-2 
from JD ~ 2452210 to 2452220 is evident in both R and 
I bands. This type of re-brightcning profile in nova evolu- 
tion is quite rare but not unique (Bonifacio et al. 2000). 
However, due to insufficient data points, we are not able to 
draw any firm conclusion about its exact behaviour. The 
evolution of (R — I) colour for the nova NMS-2 is shown 
in the lower panel of Fig. 4. It becomes bluer with time 
which is consistent with normal nova behaviour. However, 
scattering in the data due to poor flat fielding at the edge 
of frames prevents us from drawing any conclusion about 
the dust formation in the ejecta. 

The detection rate of novae in M31 is ~ 30 novae per 
year (Arp 1956, Capaccioli et al. 1989, Shatter & Irby 
2001) and 2 novae detected in our survey, keeping our field 
of view and total observing time in mind, is consistent with 
these estimates. 



5. Discussion 

The R and / magnitudes of the nova NMS-1 at the bright- 
est phase are estimated to be -7.96 and -8.02 mag respec- 
tively with a rate of decline of ~ 0.11 mag day -1 . A corre- 
lation of the rate of decline with the observed peak flux is 
in good agreement with the maximum magnitude versus 
rate of decline (MMRD) curve given by Capaccioli et al. 
(1989) for M31 novae. 

It is instructive to analyse the light curves and colour 
index, to investigate the evolution of the remnant. From 
the observed light curve of the nova NMS-1 (Fig. 2), it is 
evident that the flux in R and / bands decline in a similar 
way for the first three weeks but the late time evolution 
is qualitatively different. A difference in the profile of the 
light curves in two bands is a valuable diagnostic of the 
physical processes operating in the ejecta. We have nei- 
ther spectroscopic information nor infra-red photometry, 
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Fig. 4. R, I light curve of the nova NMS-2 during outburst 
is shown in upper panel where filled and open circles repre- 
sent R and I magnitudes respectively. The colour variation 
of the nova NMS-2 is shown in the lower panel. 



nevertheless, we might estimate the possible evolution of 
the expanding shell from the R, I light curves. 

To study the photometric behavior of the nova NMS- 
1, we divided the declining part of the light curve into 
three sections as shown in Fig. 5(al). The initial rapid 
decline in the flux follows a similar pattern in R and I 
bands and the colour remains practically unchanged; es- 
sentially this is the phase where emission from a region of 
progressively lower effective radius is received. However, 
after about three weeks, the (R — I) colour shows a grad- 
ual increase indicating that the region is cooling. After 
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J.D. 


R. 




J.D. 


R. 




t+245UUUU ) 


(mag) 


(mag) 


(+245U0U0 ) 


(mag) 


(mag) 


2198 


17.98±0.01 


17.86±0.04 


2221 




18.99±0.08 


92(12 




i 7 79 + n no 


2222 




lO.IJ ZH<J . UJ 


2203 


18.04±0.02 


17.98±0.04 


2225 


18.93±0.03 


19.11±0.06 


2206 


18.26±0.03 


18.23±0.03 


2226 


18.67±0.03 




2207 


18.40±0.02 




2229 




19.28±0.06 


2208 


18.11±0.02 




2230 


19.03±0.03 


19.41±0.07 


2213 


18.53±0.03 




2231 


18.83±0.04 


19.00±0.07 


2216 


18.43±0.03 


18.48±0.05 


2232 


19.05±0.03 


19.30±0.08 


2217 


18.42±0.02 


18.61±0.07 


2236 


19.18±0.03 




2218 


18.16±0.03 


18.20±0.04 


2237 


19.16±0.04 


19.48±0.15 
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Fig. 5. R, I light curves of the nova NMS-1 and its colour variation during the outburst is shown in the upper panels 
(al) and (a2) respectively. VJHKL data along with the colour variation of the nova Herculis 1991, taken from the 
literature, is plotted in the lower panels (bl) and (b2) respectively, to compare the behaviour of the two novae. Here 
t, to and t% indicate the time of observation, time of the maximum brightness and rate of decline respectively. 



about 40 days, the R band flux is still progressively de- 
clining similar to that observed for the nova Herculis 1991 
(Harrison & Stringfellow 1994) in the V band. The I band, 
however, probably shows a slight gradual increase in flux 
during this period; a rapid rise in the colour index of the 
nova is seen in Fig. 5(a2). A similar behaviour has been 
observed in the case of nova Herculis 1991 where an almost 
constant plateau is seen in J, H, K and L bands (see Fig. 
5(bl) and 5(b2)). Although, we have neither sufficient ob- 
servations in the later stage nor supporting spectroscopy, 
we still believe that probably there is a signature of the 
formation of neutral hydrogen. The decrease in R band 
flux but the reverse trend in I band flux is likely to be 
due to the H a absorption and re-radiation in longer wave- 



lengths. There is some reason to speculate this because: a) 
we see a nearly constant peak Balmer flux of 10 37 ergs/sec 
declining rapidly, b) moderately large velocity gradient in 
the shell is shown by Balmer lines and c) A Balmer jump 
in the spectra (cf. Downes et al. 2001) are seen. If such a 
scenario can be confirmed, it could provide a diagnostic of 
dust formation during the expansion of wind in the nova. 

For the nova NMS-2, we do not have light curve mea- 
surements in the initial phase and hence the epoch of peak 
brightness cannot be estimated. Unlike the nova NMS-1 
where we had observations for about three months, the 
nova NMS-2 could be monitored only for about 40 days 
and during this period, R and I flux show a similar trend, 
except that I declines faster. This is consistent with emis- 
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sion coming from a region of lower radius but higher tem- 
perature. Probably, the flux in both bands shows a bump 
similar to that observed in some novae. 

What does our light curves imply about the formation 
of dust in the nova wind? Novae are believed to be a source 
of dust enrichment in the interstellar medium. However, 
there are two conflicting views (cf. Bode & Evans, chap9): 
(1) Graphite or silicate grains can be formed in the cool- 
ing wind of nova depending on the chemical composition of 
white dwarf. (2) The nova environment could be hot and 
might not be the conducive of dust formation, in which 
case dust grains already present in the shell could cause 
an infra-red excess. If the expanding cool wind in the nova 
can produce grains, we should be able to see the reduction 
in temperature followed by an infra-red excess. A good 
monitoring program that tracks the temperature, lumi- 
nosity as well as possible spectroscopy will be helpful in 
settling the question. Our limited R, I observations sug- 
gest that the possible neutral hydrogen formation, when 
(R — I) increases to about 1 mag from near zero, could 
act to shield dust grains against evaporation. This might 
be the reason for the slow increase in / band flux despite 
the drop in R band flux which is contributed in part by 
regions bluewards of the Balmer alpha line. 

6. Conclusion 

We found two novae in our survey carried out during 1998 
to 2001 for an ~ 13' x 13' field belonging to the M31 disk 
population. The rate of decline of nova NMS-1 suggests 
that it was a fast nova while the decay rate for nova NMS-2 
suggests that cither it was a slow nova or we observed it at 
a very late stage after the maximum brightness phase. The 
nova NMS-1 becomes cooler with time about six weeks 
after the maximum brightness phase accompanied by a 
slow increase in / band flux, although the R band flux 
still continues to decrease. This is in contrast with the 
normal behaviour of the novae. We suspect a secondary 
bump in R, I light curves of the nova NMS-2, which is not 
unusual for normal nova evolution. 
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